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Abstract 

In this paper, we show that spatial correlation of renewable energy outputs greatly influences 
the robustness of power grids. First, we propose a new index for the spatial correlation among 
renewable energy outputs. We find that the spatial correlation of renewable energy outputs in a 
short time-scale is as weak as that caused by independent random variables and that in a long time- 
scale is as strong as that under perfect synchronization. Then, by employing the topology of the 
power grid in eastern Japan, we analyze the robustness of the power grid with spatial correlation 
of renewable energy outputs. The analysis is performed by using a realistic differential-algebraic 
equations model and the result shows that the spatial correlation of the energy resources strongly 
degrades the robustness of the power grid. Our result suggests that the spatial correlation of the 
renewable energy outputs should be taken into account when estimating the stability of power 
grids. 
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I. INTRODUCTION 


A lot of renewable energy resources are being incorporated into power networks for realiz¬ 
ing low carbon societies. More and more renewable energy resources will be introduced in all 
over the world. For example, in Japan, the amount of PV introduction in 2030 is expected 


to be about ten times larger than that in 2012 [l|. The large amount of renewable energy 
resources can cause large power fluctuations that may result in blackouts. When analyzing 
the robustness of power grids based on mathematical models, linearization schemes are not 
available in the presence of large power fluctuations, because the model behavior is far from 
a steady state. A method that enables to quantify the robustness of stable states against 
large fluctuations in differential equations was proposed in Ref. j^, where the robustness is 
measured by the volume of the basins of attraction of the stable states. This method was 
used to investigate the relation between the network topology and the robustness of power 
grids against large fluctuations j^. 

In the studies of power grids in the physics community, the Kuramoto-like phase oscillator 
model, which qualitatively corresponds to the swing equation in the electrical engineering 
community, has often been used in order to clarify various properties of power grids BQ. 


Although the model is simple and mathematically tractable, it has drawbacks that it does 
not take into account some factors which play important roles in the stability of power grids, 
particularly fluctuations in the voltage amplitude and the adjustment of reactive power in 
loads [h|. In order to study the effects of these factors on the stability of the power grids 
with relatively simple models, a new mathematical model governed by a differential-algebraic 


equations was proposec 
the swing equations [4 


121 1. In this model, the properties of generators are described using 
lOj. while the properties of loads (or substations) are determined 
with the power flow calculation. We proposed a new method to analyze the robustness 
against large fluctuations in the power grid system governed by this differential-algebraic 
equations model in Ref. jn ]. 

Among various aspects o: 


smoothing effect is notable 


renewable energy outputs intensively studied so far [l3j, the 
14J . Il5| . The smoothing effect indicates that the fluctuations 
in the total amount of renewable energy outputs are smaller than the sum of fluctuations 
in the individual renewable energy outputs because the individual fluctuations are balanced 
out. Although the smoothing effect seems to enhance the robustness of the power grid, we 
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need to take into account the spatial correlation of the weather conditions Q, Q which 
may destabilize the system locally. Therefore, it is signihcant to understand the relation 
between the spatial correlation and the robustness of power grids. 

In this study, we propose a new index that quantihes the spatial correlation of renewable 
energy resources, for characterizing the smoothing effect. In addition, we analyze the relation 
between the smoothing effect and the robustness of power grids by extending the method 
proposed in Ref. [ll|. We show that the spatial correlation, and thus the smoothing effect 
makes power grids more fragile. 

This paper is organized as follows. First, in Sec. II, we give a brief overview of the 
model of the power grids. Then, in Sec. Ill, we propose a new index for estimating the 
spatial correlation of renewable energy outputs and calculate the index value from real 
data. Section IV is devoted to clarifying the relation between the smoothing effect and the 
robustness of power grids using the real power grid topology in eastern Japan. Finally, in 
Sec. V, we conclude this paper. 


II. A DYNAMICAL MODEL OF POWER GRIDS AND A ROBUSTNESS MEA¬ 
SURE 


A. Differential-algebraic equations model 


We introduce the differential-algebraic equations model for power grids proposed by Sak- 
aguchi and Matsuo 12|]. Detailed analysis was given in Ref. 11|. 

Power grids consist of generators, loads, and end-users including buildings, factories, 
and houses. Recently, an increasing number of solar panels are connected to the end- 
users. The transmitted energy from generators to loads is equal to the difference between 
the energy consumed by end-users and the output energy of the solar panels. Therefore, 
large fluctuations in powers generated by the solar panels can cause large fluctuations of 
the transmitted energy from generators to loads. The voltages for the transmission lines 
between generators and loads are very high (500 kV or 275 kV in Japan). On the other 
hand, the voltages for the transmission lines between loads and end-users are relatively low 
(66 kV or 6.6 kV in Japan). Since the voltages are converted from high values to low values 
at substations, substations can be seen as loads from the high-voltage power grid. In this 
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study, we focus on the parts with high-voltage transmission lines, i.e., generators and loads. 
The numbers of generators, loads, and branching points are denoted by Ng, Ni, and Nf,, 
respectively, and the total number of generators and loads by 

Power grids are regarded as complex networks consisting of nodes which correspond 
to generators, loads, and branching points. We denote the voltage of node j by Vj = 
11^1 exp(i0j) for j = 1,..., A^^, where \Vj\ and 9j are the amplitude and the phase of the 
voltage, respectively. Complex numbers are useful to represent the behavior of alternate- 
current voltages. The resistance is much smaller than the reactance in the transmission 
lines, and therefore, we can consider a power network where the resistance is ne glec ted, i.e., 
a lossless network. We use the Kron’s reduction in the same way as Refs. n HQ and 
denote the reduced susceptance between node j and node k by Bj^. 

We regard the generators as phase oscillators. We denote the difference between the 
phase 9j and the phase advance by the standard frequency by = 9j — Qt. This 
standard frequency is 27r x 60 s“^ in the United States and western Japan, and 27r x 50 
s“^ in Europe and eastern Japan. The behavior of each generator is described as follows jj]: 


Nr 


Mg'cPg + \Vj BJ ^ g - (/.fc ) , 


(1) 


k=l 


where Mj is the inertial coefficient, Dj is the damping coefficient, and is the mechanical 
input energy. Equation © is called the ‘swing equation’ and derived from the energy 
conservation law. We adopt the local feedback control as follows 12l |: 


PP = -KgCj^g, 


( 2 ) 


where Kj corresponds to the inverse of the time constant. 

The behavior of each load is written as follows: 

Nr 

Fj = ^ I Vg \\Vk\ Bgk sin(0fc - , (3) 

k=l 

Nr 

Qj = \ Vg\\Vk\Bgk cos{(pk - (pj), (4) 

k=l 

where Pj and Qj are the consumed effective energy and the consumed reactive energy in 
loads, respectively. If the reactive energy is supplied in loads, Qj is negative. Equations 
(1^ and (jl]) are derived from Kirchhoff’s and Ohm’s laws. Each iteration of our numerical 
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simulation consists of two steps: First, we used the fourth-order Runge-Kutta method to 
simulate the generators’ Eq. o. We set the time step to 0.005. Second, we used Newton’s 
method for the loads’ Eqs. (E]) and (jl]). We repeated this iteration for sufficiently many 
times until a convergence is reached. We use the values of \Vj\ and at load nodes in the 
previous step to obtain (f)j of generator node j when simulating the generators’ Eq. ([T]). On 
the other hand, we use the value of (pj at generator j in the current step to obtain Pj and 
Qj of load node j when solving the loads’ Eqs. (jS]) and (jl]). 


B. Robustness evaluation 

We set the effective power and reactive power in loads stochastically. By simulating the 
generators’ Eq. ([1]) and the loads’ Eqs. ([3]) and (|4|), we judge whether the solution for Eqs. 
([3]) and (|1|) exists. The case where the solution does not exist corresponds to the power 
shortage, because this absence means that the power grid cannot realize such a power flow. 
In this case, the voltage collapse occurs. We calculate the probability that the solution exists 
and call this probability the stability rate. We evaluate the robustness based on the stability 
rate. 


III. SPATIAL CORRELATION OF RENEWABLE ENERGY OUTPUTS 


A. Smoothing effect 

The smoothing effect 1^, that the fluctuations in the total of renewable energy 
outputs are smaller than the sum of the fluctuations in the individual outputs, is a remark¬ 
able aspect of renewable energy outputs. We formulate the smoothing index in the next 
subsection. 


B. Smoothing index 

Let Tjit) and N be the power output in node j at time t and the number of nodes, 
respectively. In addition, we denote the sample standard deviation of the component of a 
row vector X by SD(X.). We quantify the strength of the smoothing effect for time window 
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size W using the following smoothing index: 


SI(W) = 


N . X 

^ SD (r,(lV + 1), r^(»' + 2),...,r,(r)) - (r,(l),rj(2 ),... ,r,(T - W))\ 

j=i V / 

. N N s 

So(5;(r,.(»y + l),r,{W + 2),.. .,r,.(r)) - ^(r,(l),rj(2 ),... ,rj{T - W))] 

Vj=i j=i / 


(5) 


where the total number of time steps is expressed as r(> W). The smoothing index dehned 
by Eq. ([S]) represents the correlation among nodes. If the correlation is zero, i.e., the 
fluctuations of noise balance out, and if probability distribution of all nodes is the same, the 
smoothing index becomes \/N. On the other hand, if rj{t) is positively correlated with each 
other, the smoothing index becomes lower than y/N. If the Pearson correlation coefficient 
(PCC) between each two nodes is 1, then the smoothing index is 1. The more the smoothing 
index is, the stronger the smoothing effect is. 


C. Smoothing index of renewable energy output data 

We calculated the smoothing index for wind speed data at 155 measurement points in 
Japan (Fig. [I]) and that for solar radiation data at 7 measurement points in the Kanto area, 
Japan (Fig. [2]) for different size of the time window. The time period for wind speed data 
and that for solar radiation data is between January 1 2010 and December 31 2012 and 
between January 1 2011 and December 31 2011, respectively. As shown in these hgures, the 
smoothing index of these data takes the value between 1 and \/#points. The dash-dotted 
line represents the case where the PCC of each two nodes is zero, i.e., where the fluctuations 
of noise are cancelled out. If we assume the independence of the fluctuations at two points, 
we obtain SI = y/ #points. The dotted line corresponding to SI = 1 represents the case 
where the PCC of each two nodes is 1. The solid line represents the smoothing index for 
the real data. In the region oi W < 12, the smoothing index decreases monotonically. 
This is because the larger the time window size is, the more similar the temporal changes 
at the individual points are, independently of their locations within Japan. This spatially 
independent trend is attributed to the nationwide weather. On the other hand, in the region 
of VP > 12, the smoothing index increases monotonically. This is because the effect of the 
daily periodicity is highest when VP = 12 and monotonically decreases when VP > 12. The 
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FIG. 1. The smoothing index of wind speed data at 155 measurement points in Japan for different 
size of time window. The horizontal axis represents the size of the time window. The step size 
is set at 10 minutes. The vertical axis represents the smoothing index. The dash-dotted line, 
indicating SI = \/l55 ~ 12.4, represents the case where the PCC of each two nodes is zero, i.e., 
where the fluctuations of noise are balanced out. The dotted line, indicating SI = 1, represents 
the case where the PCC of each two nodes is 1. The solid line represents the smoothing index for 
the real data (see the main text for details). The dashed line represents the smoothing index for 
the deviation of the data from its mean. 

output data obtained at the same time each day varies from day to day. In order to consider 
this day-by-day variability, we obtained the 24-hour data by calculating the mean of the 
output data recorded at the same time each day over the days of observation. The deviation 
of the output data from its mean is represented by the difference between the output data 
of each day and the 24-hour data. The dashed line represents the smoothing index of the 
deviation of the output data. The difference between the solid line and the dashed line 
represents the contribution of the daily periodicity. The deviation of the output data from 
its mean is more representative of the smoothing effect. 

In summary, the smoothing index of the renewable energy outputs decreases in time. This 
implies that the renewable energy outputs of different points fluctuate independently in a 
short time scale, and the dependence is stronger if we consider a longer time window. This 
result suggests the necessity to study stability of power grids taking into account the spatial 
correlation of the renewable energy outputs. In the next section, we introduce a spatial 
correlation of the effective power of different sites in the model and study dependence of the 
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FIG. 2. The smoothing index for solar radiation data at 7 measurement points in Japan for different 


size of the time window. The horizontal axis represents the size of the time window. The step size 


is 1 minute. The vertical axis represents the smoothing index. The dash-dotted line, indicating 
SI = \/7, represents the case where the PCC of each two nodes is zero, i.e., where the fluctuations 


of noise are balanced out. The dotted line, indicating SI = 1, represents the case where the PCC 


of each two nodes is 1. The solid line represents the smoothing index for the real data. The 


dashed line represents the smoothing index for the deviation of the data from its mean. In the 


solar radiation data, there exist corrupted data due to the breakdown or checkback of measurement 


equipments. We excluded such data when calculating the mean and the smoothing index. 


stability of the power grid on the spatial correlation. 


IV. IMPACT OF SMOOTHING EFFECT ON ROBUSTNESS OF THE POWER 
GRIDS 


We simulated the differential-algebraic equations model for the power grid in eastern 
Japan by taking into account the smoothing effect found in the previous section. 













FIG. 3. (Color online) The topology of the high voltage power grid in eastern Japan that we 
derived. The red circles represent generators, the green circles represent loads, and the small gray 
circles represent branching points. 


A. Eastern Japan power grid 


Figure E] shows the topology of the high voltage power grid in eastern Japan. It was 


derived from the reports released 


Electric Power Company Inc. 


by Tokyo Electric Power Company Inc. 


20 1 and Tohoku 


2l| . The hgure was produced with the graph visualization 


software Gephi [22|. This power grid consists of Ng = 53, Ni = 126, and Nb = 69 nodes. 
The number of nodes corresponding to generators and loads is Nr = Ng + Ni = 179 and 
the total number of nodes in the power grid is N = 248. In addition, we set the parameter 
values to Mj = 1, Dj = 1, Kj = 30, PJ”(0) = 1, and Qj = 0. The admittance matrix for the 
links was set as follows jll| : 


Bjk — 


10, if node j and node k are connected. 


0 , 


otherwise. 


( 6 ) 
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B. Results 



FIG. 4. The relation between the variance of fluctuations and the stability rate. The horizontal axis 
represents the variance of the fluctuations in each node. The vertical axis represents the stability 
rate. The open circles with the solid line show the result when the PCC of each two nodes is 0. 
The asterisks with the dashed line correspond to the result when the PCC of each two nodes is set 
to 0.5. 


We give a fluctuation of consumed effective power in each load, i.e., Pj. The mean of 
the effective power Pj is set to 0.45. In addition, we vary the variance of Pj from 0 to 
0.15. For simplicity, we assume that the mean and the variance of the fluctuations are 
identical for all nodes and the PCC of each two nodes is the same. Further, we assume 
that the autocorrelation of the fluctuations in each node is 0. Figure H] shows the relation 
between the variance of fluctuations and the stability rate. The open circles with the solid 
line indicate the result when the PCC of each two nodes is 0. The asterisks with the dashed 
line indicate the result when the PCCs are set at 0.5. Clearly, the power grids become more 
unstable when the fluctuations are more correlated. 

We further hxed the variance of fluctuations in each node and varied the PCC of those 
between each two nodes. Figure [5] shows the relation between the PCC of each two nodes 
and the smoothing index. The vertical axis represents the smoothing index dehned in Eq. 
(]5|). The time window size W is set to be equal to one step of the numerical simulation. The 
horizontal axis represents the PCC of fluctuations between each two nodes. The larger the 
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FIG. 5. The relation between the PCC of fluctuations between each two nodes and the smoothing 
index. The vertical axis represents the smoothing index defined in Eq. ([5]). The time window size 
W is equal to one step of the numerical simulation. The horizontal axis represents the PCC of 
each two nodes. The open circles with solid lines show the experimental result. The asterisks with 
dashed lines show the analytical solution, which is derived in Appendix. It is in good agreement 
with the experiment result. 

PCC is, the smaller the smoothing index is. When the output of each node is independent, 
the smoothing index is \/W = \/120. In fact, the smoothing index is near \/M when 
the PCC is 0. Figure [6] shows the relation between the PCCs and the stability rate. The 
horizontal axis represents the PCC of each two nodes. The vertical axis represents the 
stability rate. The larger the PCC is, the less stable the power grid is. Figure [7] shows the 
relation between the smoothing index and the stability rate. The horizontal axis represents 
the smoothing index dehned in Eq. ([S]). The time window size W was set to be one step 
of the numerical simulation. The vertical axis represents the stability rate. The larger 
the smoothing index is, the more stable the power grid is. These results mean that the 
smoothing effect can help stabilize the power grid. 


V. CONCLUSIONS 

In this paper, we have proposed a new index, called the smoothing index, which quantihes 
the spatial correlation among renewable energy outputs. We have analyzed the relation 
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FIG. 6. The relation between the PCC of fluctuations and the stability rate. The horizontal axis 
represents the PCC of each two nodes. The vertical axis represents the stability rate. 



FIG. 7. The relation between the smoothing index and the stability rate. The horizontal axis 
represents the smoothing index defined in Eq. ([5]). The time window size W was set to be equal 
to one step of the numerical simulation. The vertical axis represents the stability rate. 


between the smoothing effect and the robustness of power grids by using the mathematical 


model governed by differential-algebraic equations proposed in Ref. 


12| and the real topology 


of the power grid in eastern Japan. We have clarified that the smoothing effect facilitates 
the stability of the power grid. In other words, the spatial correlation of the renewable 
energy outputs degrades the robustness of the power grid, which is an important factor for 
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studying the introduction of the renewable energy. 
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APPENDIX: DERIVATION OF THE RELATION BETWEEN THE SMOOTH¬ 
ING INDEX AND THE PCC 

For simplicity, we assume in the simulation that the mean and the variance of fluctuations 
are identical for all the nodes, the PCC between any pair of two nodes is the same, and the 
autocorrelation of each node is 0. In the approximate derivation of the relation between the 
smoothing index and the PCC, we assume that the sample standard deviation is equal to 
the standard deviation and the correlation between rj(t + 1) and rk{t) is 0 for simplicity. 
By setting the variance V{rj{t)) to and the correlation of each node to c, the smoothing 
index is calculated as follows: 

MEf=i(c(^ + i)-cW)) 

EjLi + 1) - rj(Jj) 

NV(V(n{t + l)) + V(,r,m 

nV^ 

~ ^2Na^ + 2N{N - l)ca‘^ 

= , ^ . (7) 

yi + (Af-l)c 

Therefore, *5'/(l) = \fN if c = 0 and SI{1) = 1 if c = 1. The analytic solution of the relation 
between the smoothing index «S'/(1) and the correlation is depicted for real data in Fig. [5l 


13 














It is in good agreement with the experimental result. 
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